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ABSTRACT 

In reionized regions of the Universe, gas can only collapse to form stars in dark matter 
(DM) haloes which grow to be sufficiently massive. If star formation is prevented in 
the minihalo progenitors of such DM haloes at redshifts z > 20, then these haloes will 
not be self-enriched with metals and so may host Population (Pop) III star formation. 
We estimate an upper limit for the abundance of Pop III star clusters which thus 
form in the reionized Universe, as a function of redshift. Depending on the minimum 
DM halo mass for star formation, between of the order of one and of the order of a 
thousand Pop III star clusters per square degree may be observable at 2 <^ z <, 7. 
Thus, there may be a sufficient number density of Pop III star clusters for detection 
in surveys such as the Deep- Wide Survey (DWS) to be conducted by the James Webb 
Space Telescope (JWST). We predict that Pop III clusters formed after reionization 
are most likely to be found at z i; 3 and within ~ 40 arcsec (~ 1 Mpc comoving) of 
DM haloes with masses of ~ 10 11 M Q , the descendants of the haloes at z ~ 20 which 
host the first galaxies that begin reionization. However, if star formation is inefficient 
in the haloes hosting Pop III clusters due to the photoionizing background radiation, 
these clusters may not be bright enough for detection by the Near-Infrared Camera 
which will conduct the DWS. Nonetheless, if the stellar initial mass function (IMF) 
is top-heavy the clusters may have sufficiently high luminosities in both Lya and 
He II A1640 to be detected and for constraints to be placed on the Pop III IMF. While 
a small fraction of DM haloes with masses as high as ~ 10 9 M at redshifts z 4 
are not enriched due to star formation in their progenitors, external metal enrichment 
due to galactic winds is likely to preclude Pop III star formation in a large fraction of 
otherwise unenriched haloes, perhaps even preventing star formation in pristine haloes 
altogether after reionization is complete at z ~ 6. 

Key words: cosmology: theory - early Universe - galaxies: formation - high-redshift 
- haloes - intergalactic medium 



1 INTRODUCTION 

How long did the epoch of primordial star formation last? 
This is an important and timely question, as a primary goal 
of next generation telescopes, such as the James Webb Space 
Telescope (JWST; e.g. Gardner et al. 2006; Windhorst et 
al. 2006) and the European Extremely Large Telescope^ (E- 
ELT), is to detect primordial galaxies and star clusters in 
the early Universe. While these telescopes may allow to de- 
tect individual supernovae (SNe) exploding within the first 
Population (Pop) III star clusters and galaxies (see e.g. Wise 
& Abel 2005; Haiman 2008), the emission from the stellar 
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populations within the first dwarf galaxies at z J; 10 is likely 
too dim for even the JWST to detect (e.g. Ricotti et al. 2008; 
Johnson et al. 2009). Instead, only very rare, massive pri- 
mordial galaxies at such high redshifts, or dwarf galaxies and 
star clusters forming at lower redshift may be bright enough 
for detection (see e.g. Ciardi & Ferrara 2001; Scannapieco et 
al. 2003). Therefore, whether or not Pop III star clusters can 
be found by next generation telescopes critically depends on 
the amount of primordial star formation that takes place at 
relatively low redshift. 

The hierarchical nature of structure formation dictates 
that more massive dark matter (DM) haloes are built up 
from smaller ones, and that more massive haloes generally 
have progenitors in which star formation may take place at 
earlier times. In order for a Pop III star cluster to form in a 
DM halo, star formation in its progenitor haloes may have to 
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be completely suppressed, lest SNe enrich the primordial gas 
with heavy elements (but see e.g. Jimenez & Haiman 2006; 
Pan & Scalo 2007; Dijkstra & Wyithe 2007). Exposure to 
a persistant photoionizing background, such as takes place 
in reionized regions of the Universe at high redshifts, can 
suppress star formation in DM haloes with masses substan- 
tially higher than those in which the first stars (e.g. Abel 
et al. 2002; Bromm et al. 2002; Yoshida et al. 2003) and 
galaxies (e.g. Greif et al. 2008; Bromm et al. 2009) form, 
completely preventing the infall of gas into haloes with cir- 
cular velocities up to ~ 30 km s~\ corresponding to virial 
masses of £ 3 x 10 8 M at z <, 10 (e.g. Efstathiou 1992; 
Thoul & Weinberg 1996; see also Kitayama & Ikeuchi 2000; 
Gnedin 2000; Ciardi & Ferrara 2005). If star formation in the 
progenitors of these more massive haloes can be suppressed 
by the photoionizing background during reionization, then 
these haloes can remain unenriched until the primordial gas 
collapses into them to form clusters of Pop III stars at rela- 
tively low redshifts, e.g. at z Si 6. 

The formation of Pop III stars at such low redshifts 
has been a possibility recently studied through the use of 
large-scale cosmological simulations. In their high-resolution 
simulations of cosmic reionization, Shin et al. (2008) (also 
Trac & Cen 2007) found that haloes with masses ~ 10 8 
Mq may host Pop III star formation at z ^ 6, but their 
simulation did not resolve the collapse of minihaloes, i.e. 
haloes with virial temperature T v \ r 10 4 K, which could 
have led to the enrichment of the gas at higher redshifts (see 
also Schneider et al. 2006; Tornatore et al. 2007). In another 
recent study, Trenti et al. (2009) model the enrichment of 
haloes with T v i r <; 10 4 K by their progenitors, and find that 
significant Pop III star formation may take place at z Si 6, as 
well. However, their study did not explicitly account for the 
suppression of star formation by photoionization, although, 
especially at redshifts z <; 10 (Dijkstra et al. 2004), it is 
possible that star formation takes place in haloes with T V1I 
~ 10 4 K even in the presence of a photoionizing background. 

In related analytical work, Wyithe & Cen (2007) esti- 
mated the contribution of Pop III stars to the reionization of 
the Universe, arguing that primordial star formation in neu- 
tral minihaloes may be suppressed until these haloes grow 
to be sufficiently massive, perhaps allowing for the forma- 
tion of a significant fraction of Pop III stars up to the end 
of reionization at z ~ 6. In this work, it was assumed that 
star formation in minihaloes was prevented by the photodis- 
sociation of molecules at very high redshifts, e.g. z <; 20 
(Haiman et al. 1997). However, more recent work has shown 
that molecular cooling can operate in minihaloes despite the 
build-up of the photodissociating background radiation field 
(e.g. Wise & Abel 2007; O'Shea & Norman 2008; Johnson et 
al. 2008; Trenti & Stiavelli 2009), providing motivation to 
relax this assumption and consider alternate scenarios for 
the formation of Pop III stars at low redshifts. 

In the present work, we account for star formation in 
minihaloes at z <; 20, as well as for the suppression of star 
formation in DM haloes during reionization, with the aim of 
evaluating the abundance and detectability of Pop III star 
clusters in reionized regions of the Universe. In Section 2, we 
outline the requirements for Pop III star cluster formation 
in the post-reionization Universe. In Section 3, we calcu- 
late an upper limit for the abundance of such clusters and 
estimate the effect that large-scale metal-enriched galactic 



winds have in preventing their formation, while in Section 4 
we discuss the prospects for their detection in future surveys. 
We conclude with a discussion of our results in Section 5. 
For our calculations, we assume the latest cosmological pa- 
rameters reported by the Wilkinson Microwave Anisotropy 
Probe (WMAP) collaboration (Komatsu et al. 2009). 



2 CONDITIONS FOR POPULATION III STAR 
FORMATION AFTER REIONIZATION 

The first stars are expected to have formed in DM mini- 
haloes with virial temperature ~ 10 3 K at redshifts z <; 20 
(see e.g. Bromm & Larson 2004; O'Shea & Norman 2007). 
These stars were likely very massive, with a fraction of them 
ending their lives as powerful supernovae which enriched the 
intergalactic medium (IGM) with the first heavy chemical 
elements (e.g. Bromm et al. 2003; Kitayama et al. 2005; 
Greif et al. 2007; but see also Whalen et al. 2008a). For a 
Pop III star cluster to form in a given DM halo, the forma- 
tion of stars which explode as SNe and enrich the gas may 
have to be completely prevented in its minihalo progenitors. 
We take this as a first condition for Pop III star formation 
in the reionized Universe. For simplicity, we make the con- 
servative assumption that no Pop III stars can form in the 
progenitors of a given halo in order for it to remain pristine, 
although in principle primordial stars which collapse to form 
black holes, instead of exploding as SNe, may not enrich the 
gas in their host haloes (see Fryer et al. 2001; Heger et al. 
2003; Schneider et al. 2006). We discuss the related issue 
of the metal enrichment of haloes by large-scale winds from 
neighboring galaxies in Section 3.2. 

A molecule-dissociating, so-called Lyman- Werner (LW), 
background radiation field can delay star formation in dark 
matter (DM) minihaloes, but the gas eventually collapses 
when the DM halo becomes massive enough, of the order of 
10 7 M for very high LW flux (Shang et al. 2009; see also 
Wise & Abel 2007a; O'Shea & Norman 2008; Mesinger et 
al. 2006, 2009). Furthermore, as we show below, star forma- 
tion in minihalo progenitors at z <; 19, on average, must be 
suppressed for haloes forming in the reionized Universe to 
avoid self-enrichment by these progenitors. Recent modeling 
of the build-up of the LW background suggests that the LW 
flux may remain too low to prevent star formation in all but 
the rarest minihalos at z ^ 15 (Johnson et al. 2008; Dijkstra 
et al. 2008; Ahn et al. 2009; Trenti & Stiavelli 2009). There- 
fore, we expect that LW feedback alone is not sufficient to 
prevent the collapse of gas into the progenitors of DM haloes 
with masses J> 10 8 M , the mass range for haloes in which 
star formation can take after reionization is complete at z 
~ 6. 

A more potent form of radiative feedback at high red- 
shifts is the photoionizing background radiation that builds 
up during reionization. In regions which are reionized early, 
at e.g. z <; 15, the temperature of the photoheated gas is too 
high, at <; 10 4 K, for minihaloes to retain gas, and star for- 
mation can be suppressed. We note that for star formation 
to be suppressed by the photoionizing background during 
reionization, the minihalo must be subjected to the radia- 
tion before it reaches advanced stages of collapse, as an ion- 
ization front passing through a dense minihalo may act to 
induce star formation, not to prevent it (e.g. Ahn & Shapiro 
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Figure 1. The maximum mass M max (z r cion ) of DM haloes which have, on average, no minihalo progenitors massive enough to host star 
formation before a redshift z r cion> as a function of redshift z. Also shown are two possible minimum halo masses, M m ; n 20 an d A^min,30i 
required for photoionized gas to collapse and form stars. For minimum halo masses of M m j n 20 an< i ^min,30i on ly haloes in which star 
formation is suppressed by the reionization of the gas at z r cion > 19 and > 23, respectively, can remain pristine until the reionized gas 
collapses to form stars. These star-forming haloes may host Pop III star clusters luminous enough to be detected in deep surveys by the 
JWST (see Section 4). 



2007; Whalen et al. 2008b; Susa et al. 2009). Therefore, we 
make the simple, conservative assumption that the gas col- 
lapsing into minihalos with masses higher than M m i n ji m io = 
10 6 M ((1 + z)/10)~ 3 '" 2 , corresponding to virial tempera- 
ture T vir ~ 2 x 10 3 K (see e.g. Yoshida et al. 2003), must be 
reionized before the virialization of the minihalo, in order for 
the collapse of the gas and star formation to be prevented. 
This assumption could be relaxed slightly if we account for 
the photoevaporation of the gas inside minihaloes that is 
already somewhat collapsed, but we do not expect this to 
dramatically affect our results (see e.g. Shapiro et al. 2004; 
Iliev et al. 2005; Wyithe & Cen 2007). 

A second condition for Pop III star formation in unen- 
riched DM haloes after reionization is that the DM haloes 
must be massive enough to allow the photoheated gas to cool 
and collapse. Numerous studies have found that DM haloes 
must have circular velocities <; 30 km s~ in order for pho- 
toionized gas to collapse into them (e.g. Thoul & Weinberg 
1996; Quinn et al. 1996; Kepner et al. 1997; Kitayama et al. 
2001). This limit may be lower at redshifts z <; 10, due to 
the lower average intensity of the ionizing background and 
to the higher densities to which baryons may have already 
collapsed when the ionizing radiation turns on (Kitayama 
& Ikeuchi 2000; Dijkstra et al. 2004; see also Chiba & Nath 
1994). However, there is no concensus on what this limit- 
ing mass is in cosmological haloes (Wyithe & Cen 2007), 
and many previous calculations have been done assuming 
spherical symmetry (but see the cosmological simulations of 
Hoeft et al. 2006), which may underestimate the minimum 
halo mass, as asphericity can allow ionizing photons to pene- 
trate deeper into haloes and reduce the effect of self-shielding 
(see Haiman 2008). However, this is not necessarily so clear 



either, as dense filaments from which gas can accrete into 
cosmological haloes may also be self-shielding, in which case 
calculations in sperical symmetry would lead to an overes- 
timate of the minimum mass. To account for some spread 
in the minimum mass required for gas to collapse into DM 
haloes, we consider two possibilities for this minimum mass, 
corresponding to minimum circular velocities of 20 and 30 
km s _1 , given below: 

M min , 20 = 9 x 10 7 M (1) 

M min , 30 = 3 x 10 s M (-^) ■ (2) 

Using the extended Press- Schechter formalism (Lacey 
& Cole 1993), we calculate, as a function of redshift, the 
maximum mass M max (2: re ion) of virialized DM haloes which 
have, on average, no minihalo progenitors with masses J> 
Afminihaio prior to redshift £ re ion; more precisely, we calculate 
the maximum mass of haloes which have eW/dlnMmmhaio < 
1, where N is the number of progenitor haloes at z rc ion (see 
e.g. Bromm & Clarke 2002). This nomenclature is chosen 
to denote the redshift by which the region in which a halo 
forms must be reionized in order to prevent star formation 
in these minihalo progenitors. 

We compare this maximum mass, for five different 
choices of z rc i on , to the minimum mass for star formation 
in Figure 1. Only in haloes with masses below the maxi- 
mum halo mass M max (z rc i on ) and above the minimum mass, 
-Mmi nj 20 or M m i nj 3o, at a given redshift z can the reionized 
gas remain pristine and unenriched until it collapses to form 
stars. Accordingly, in regions reionized at z rc i on ~ 25, pri- 
mordial star formation may take place at z 12.5, on av- 
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Figure 2. A schematic illustration of the scenario for Pop III star cluster formation in reionized regions of the Universe described in 
Sections 2 and 3. Star formation in minihaloes which are massive enough to host Pop III star formation in the absence of radiative 
feedback (~ 10 6 Mq) is suppressed at redshifts z < 20 by photoheating of the primordial gas within the early H II region (dark gray; 
with a typical comoving radius of < 1 Mpc) created by a dwarf galaxy formed in a halo with a virial temperature > 10 4 K (light blue; 
corresponding to a mass of ~ 10 8 Mq at z ~ 20). The minihalo at z ~ 20 grows to a mass of M m ; n , the minimum required for star 
formation in a photoheated gas, at z < 6 (see Section 2), at which point the descendant of the first galaxy at z ~ 20 would typically be 
hosted in a ~ 10 11 Mq halo (see Section 4). This galaxy is likely to drive a metal-enriched wind (light gray) into the IGM at some velocity 
^wind ( see Section 3.2). If this wind does not reach the halo with mass M m ; n by the time stars are formed, or if the metals transported by 
the wind are not mixed into the halo, then the halo is assumed to host a Pop III star cluster. Note that while the approximate size of the 
H II region and the size of the metal-enriched bubble relative to it are roughly to scale (although shown as spherical here for simplicity), 
for clarity the galaxy, cluster, and halo sizes are not shown to scale. 



erage, if the minimum halo mass is M m i nj 2o- If, instead, the 
minimum halo mass is M m j n ,3o then primordial star forma- 
tion can only take place at z ^ 8. Also, the most massive 
haloes that can host metal-free star formation form at z 4 
in regions that are reionized at z rc ion = 19 and 23, for a min- 
imum halo mass of M m i nj 20 or M m i nj 3o, respectively. Indeed, 
due to the increase in the minimum mass for star formation 
with decreasing redshift, it is only in regions reionized at z 
<; 19 in which Pop III star formation can take place in reion- 
ized regions of the Universe, unless the minimum halo mass 
for star formation is considerably smaller than M m i n ,2Q. 

To summarize, two minimum requirements for the for- 
mation of a Pop III star cluster in an unenriched DM halo 
after reionization, i.e. at redshifts z 6, are the following 
(see also Fig. 2): 

(1) The halo must not have any progenitors which host star 
formation. This implies a maximum halo mass M max (z rc i on ), 
above which star formation occurs in progenitor minihaloes 
that form prior to the redshift z rc i on at which the gas is 
reionized. After z re ion, photoheating of the gas prevents star 
formation, and the gas remains unenriched. 

(2) The halo must be massive enough for photoheated gas 
to collapse into it and form stars. This implies a minimum 
halo mass M m i n , below which stars are unable to form. 

We thus take it that only haloes with masses <; M m i n 
and < M max can host the formation of primordial stars in 



reionized regions of the Universe. A third requirement is that 
the haloes are not enriched by galactic winds emanating 
from neighboring star-forming haloes, as shown schemati- 
cally in Figure 2. Such external chemical enrichment will 
limit the abundance of unenriched haloes; this is discussed 
further in the next Section. 



3 THE ABUNDANCE OF POP III STAR 
CLUSTERS AFTER REIONIZATION 

3.1 An upper limit 

Now that we have the redshift and mass distribution of DM 
haloes which may potentially host Pop III star formation 
in reionized regions of the Universe, we would like to esti- 
mate the abundance of such haloes as a function of redshift, 
thereby finding an upper limit for the abundance of Pop III 
star clusters. As shown in Section 2, these haloes can only 
form in the relatively rare regions of the Universe that were 
reionized at very early times, i.e. at z rc ion 19, meaning 
that these haloes are relatively rare objects. While the early 
reionization history is only poorly constrained by observa- 
tional data, there have recently been conducted a number of 
large-scale high-resolution simulations of this process which 
yield useful results for our purposes (e.g. Iliev et al. 2007; 
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Shin et al. 2008). Shin et al. (2008) report their result for 
the volume fraction / r cion,voi(z r cion) of the Universe that is 
reionized, as a function of redshift, up to z rc i on ~ 22. For our 
calculation, we use the following formula, which provides a 
satisfactory fit to their results over the range of redshifts 
Zroion that we consider (i.e. 19 Si z re ion S> 25): 



log(/, 



rcion f vol 



/ \ \ t , / ^ 1 ^reion \ 
(Zrcion)) = -4 + " J 



(3) 



The rate of change of the collapse fraction fni(z rc ion, z) 
of mass in haloes which may host Pop III star formation at 
redshift z, if forming in regions that are reionized at z re ion, 
is given by the rate at which mass becomes incorporated in 
haloes with the minimum mass for star formation, as follows: 
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if M min < M max ( ^reion) ; 
otherwise. 



(4) 



Here F(M,z) is the collapse fraction in DM haloes with 
masses > M at redshift z, as given by the Sheth-Tormen 
formalism (Sheth & Tormen 1999). M m i n is the minimum 
mass for star formation, given either by equation (1) or (2). 

As the Sheth-Tormen formalism gives the fraction of 
mass that is collapsed, it is more appropriate to use the 
fraction of mass / rc ion,mass that is reionized than the volume- 
filling fraction / ro ion,voi for our calculation. We shall assume 
that the mass fraction of the Universe that is reionized is a 
factor of ten higher than the volume- filling fraction / ro ion,voi 
given in equation (3) at 19 Si z <> 25, due to the relatively 
large overdensities of the rare regions that are reionized so 
early. This is roughly consistent with the results of Shin 
et al. (2008), who report that at z ~ 20 the volume-filling 
fraction of H n regions is of the order of 10 -4 while the 
fraction of the gas that is ionized is of the order of 10~ 3 , 
implying an average overdensity of the H n regions of the 
order often (see also Furlanetto et al. 2004). Accounting for 
this fraction / rc ion,mass(z ro ion) of the baryonic mass in the 
Universe that was reionized before z rc ion, we integrate over 
Zrcion to find the total number density of haloes virializing at 
redshift z which are unenriched by previous star formation. 
Assuming that a single Pop III star cluster with a lifetime 
of tdustcr forms in every unenriched halo, we thus estimate 
the comoving number density of Pop III clusters as 
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where p cr i t is the critical density of the Universe, assuming 
a Hubble constant Ho = 70.5 km s _1 Mpc -1 (Komatsu et 
al. 2009), and Qdm is the fraction of this density in dark 
matter. The integral is taken from z re i n = 19 to 25, because 
unenriched haloes with masses <; M m i n only form in regions 
reionized before Zrcion = 19, as found in the previous Section, 
and because before z rc ion = 25 such a small fraction of the 
Universe is reionized that little contribution would be made 



to the integral for higher z re i n. Finally, we divide by the 
minimum halo mass to find the number density of haloes. 

We plot two curves for the number density nm of 
Pop III star clusters in the top panel of Figure 3, one for each 
of the choices of M m i n , given by equation (1) and (2). The 
more massive haloes with mass M m i nj 3o are much rarer than 
those with mass M m i n ,20 for two reasons. The first is that 
the overall abundance of DM haloes generally decreases with 
halo mass. The second reason is that more massive haloes 
have minihalo progenitors at higher redshift, meaning that 
on average they form only in the regions in which star for- 
mation is suppressed by reionization at higher z rc i n, and 
which are correspondingly more rare. For reasons discussed 
in Dijkstra et al. (2004), the lower minimum mass may be 
more accurate at higher redshift, z Sj 10, while at redshifts 
z Si 3 the higher minimum mass is likely more accurate. If 
this is so, then the comoving number density of Pop III star 
clusters may remain below ~ 10 4 Gpc~ 3 at all redshifts. 

This would be consistent with the upper limit on the 
abundance of Pop III clusters at z ~ 4 found by Nagao et 
al. (2008), shown in the top panel of Fig. 3. These authors 
report the non-detection of strong emitters of both Lya and 
He n A1640 within a comoving volume of 4.03 x 10 5 Mpc 3 
in the Subaru Deep Field at 3.93 Ss z Si 4.01 and 4.57 Ss z 
<, 4.65. This corresponds to an upper limit on the number 
density of Pop III clusters of tim ~ 2500 Gpc -3 , between 
the upper limits plotted in Fig. 3 for minimum halo masses 
of M m i ni 2o and M m i n ,3o. While this upper limit suggests that 
either the minimum halo mass for star formation is higher 
than M m i nj 2o at low redshifts or that metal enrichment by 
galactic winds is relatively rapid (see Section 3.2), we note 
that the limiting flux of the survey would have allowed to 
detect only the brightest clusters (see Section 4). 

In the bottom panel of Fig. 3, we plot the number den- 
sity of Pop III star clusters per degree observed on the sky 
per redshift interval Az. We also plot the inverse of the 100 
arcmin 2 area of the planned Deep- Wide Survey (DWS) to be 
conducted by the JWST, in which Pop III star clusters may 
be detected (Gardner et al. 2006); this corresponds to the 
number density of clusters that is required for there to be 
at least one cluster within this survey area, on average, per 
unit redshift. The number of clusters that the JWST may 
detect is clearly very sensitive to the minimum halo mass for 
star formation; most optimistically, for M m i ni 2o, there may 
be up to ~ 100 clusters out to z ~ 6 per DWS area. How- 
ever, for a minimum halo mass of M m i n ,3o the total number 
density is too low for the detection of even a single cluster 
in such a survey. A larger survey area of, for example, 17 
x 17 arcmin 2 for the E-ELT (ESO ELT Science Working 

Group t) would, of course, provide a better chance to detect 
Pop III clusters. 

As shown in equation (5), the number density of Pop III 
clusters is linearly proportional to the lifetime truster of an 
individual cluster. Thus one impediment to the detection 
of Pop III clusters is their brief lifetime. Because a halo 
becomes enriched with metals once the first SNe explode, 
the lifetime of a purely primordial cluster is in part deter- 
mined by the lifetime of the massive stars which explode as 



t http://www.eso.org/sci/facilities/eelt/docs/ELT-SWG-apr30- 
l.pdf 
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Figure 3. Upper limits for the number density of Pop III star clusters nin formed from reionized gas, as given by equation (5), for 
our two choices for the minimum DM halo mass for star formation given by equations (1) and (2). Here nni is normalized to a cluster 
lifetime t c i us tcr = 2 X 10 s yr, the expected lifetime of the most massive Pop III stars. Top panel: the number of clusters per comoving 
Gpc — 3 . Also shown is the upper limit to the number density of bright Pop III star clusters at z ~ 4 reported by Nagao et al. (2008) (see 
text in Sections 3 and 4). Bottom panel: the number of clusters per redshift interval Az per square degree on the sky. The dotted line 
shows the number density of clusters at which there is one cluster per redshift interval, on average, within the 100 arcmin 2 area of the 
JWST Deep- Wide Survey. This survey could reveal at most ~ 100 clusters out to z ~ 6, but this depends sensitively on the minimum 
DM halo mass for star formation in a reionized gas; if the minimum circular velocity for such haloes is 30 km s — 1 instead of 20 km a -1 , 
then such a survey would likely turn up no Pop III clusters. 



SNe, which is, at most, of the order of 10 7 yrs (Schaerer 
2002). However, the timescale for the formation of second- 
generation, Pop II stars from the enriched gas is also de- 
pendent on the timescale for the mixing of metal-enriched 
SNe ejecta with the primordial gas, which could be J> 10 8 
yr (see e.g. de Avillez & Mac Low 2002). In this case, the 
timescale for the formation of enriched stars would be simi- 
larly long; accordingly, the abundance of DM haloes hosting 
only Pop III stars would be /& 50 times higher than in the 
fiducial case of Cluster — 2 x 10 6 yr shown in Fig. 3. We 
emphasize, though, that it is the most massive stars which 
are the most luminous and hence it is during the lifetime 
of these stars, a few times 10 6 yrs, that a cluster is most 
easily detectable (see Section 4) , meaning that the brightest 
clusters are also likely to be the rarest. 

While in Figure 3 we have only plotted the Pop III 
cluster abundances for one value of the minimum halo mass 
M minih alo = 10 6 M ((1 + z)/10) _3/2 required for star for- 



mation in minihalos, this minimum halo mass is an uncer- 
tain quantity and may be higher if a persistant LW back- 
ground radiation field is present (e.g. O'Shea & Norman 
2008; Trenti & Stiavelli 2009). If M miniha i is indeed higher, 
then the abundance of Pop III clusters would be higher as 
well, since these potentially star-forming progenitor mini- 
haloes would form at lower redshifts when /reion.mass is con- 
siderably higher. For a minimum minihalo mass an order of 
magnitude higher, M mini haio = 10 7 M Q ((1 + z)/10)" 3/2 , we 
find that the upper limits on the abundance of Pop III star 
clusters may be roughly an order of magnitude higher than 
those shown in Fig. 3. 

3.2 External metal enrichment 

We have so far considered the formation of star clusters in 
DM haloes which are unenriched by star formation in their 
progenitors. However, metal enrichment can also take place 
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Figure 4. The number density rani of Pop III star clusters, taking into account enrichment of the IGM by galactic winds. As in Fig. 
2, the solid lines correspond to a minimum halo circular velocity of 20 km s — 1 , while the dashed lines correspond to a minimum of 30 
km s — 1 . For each series of lines, the top (black) line is our model which neglects external metal enrichment (same as in Fig. 3), the green 
line corresponds to a wind velocity f w i n< j = 30 km s — , the red line to f w i n( j = 60 km s , and the blue line to t> w i n d = 100 km s~ 1 . For 
wind velocities of ~ 60 km s" 1 , the volume- filling fraction of metals in the early H II regions in which Pop III clusters may later form 
reaches unity at z ~ 3, thereby shutting off Pop III star formation altogether. For wind velocities higher than 100 km s , this occurs 
at z > 5. 



via galactic winds launched by neighboring star-forming 
galaxies (e.g. Nath & Trentham 1997; Mac Low & Ferrara 
1999; Cen & Bryan 2001; Madau et al. 2001; Mori et al. 2002; 
Scannapieco et al. 2003; Bertone et al. 2005; Tornatore et 
al. 2009). As such external metal enrichment was ignored in 
the last Section, our results so far for the number density of 
Pop III clusters are only upper limits. To estimate the extent 
to which galactic winds enrich the IGM and thereby prevent 
Pop III star formation, we compute the volume-filling frac- 
tion / w i n d of the IGM that is reionized at z 19 which these 
winds sweep up, as a function of redshift. 

As we are interested in the metal enrichment of reion- 
ized regions, we assume that winds only originate from DM 
haloes with masses greater than the minimum mass for star 
formation in reionized regions. While haloes with masses 
~ Af,nin will be the much more abundant than more mas- 
sive haloes (see e.g. Reed et al. 2007), the winds emanating 
from the galaxies which form in these haloes are likely to be 
rather weak, as the star formation efficiency in these haloes 
is low (e.g. Furlanetto & Loeb 2003; Samui et al. 2008), due 
to the haloes retaining only a rather small fraction (~ 10 
- 20 percent) of their gas due to photoheating (e.g. Thoul 
& Wienberg 1996; see also Shen et al. 2009). However, the 
regions of interest to us are in close proximity to the descen- 
dants of the first galaxies which reionized their surroundings 
at Zreion <; 20. As these first galaxies formed in haloes with 
masses of £ 5 x 10 7 M at z £ 20 (Trac & Cen 2007), their 
descendants, on average, are rare haloes with masses <; 10 10 
Mq at z 10, as found using the extended Press-Schechter 
formalism. These galaxies likely began enriching the IGM 
already at z ~ 20, and because star formation is likely to 
proceed with a high efficiency within such massive haloes, we 



expect that these galaxies launch the metal-enriched winds 
which are principally responsible for enriching the otherwise 
pristine haloes we are considering. Thus, under this assump- 
tion, as described schematically in Figure 2, it is the same 
galaxies which reionize the primordial gas at early times, 
thereby suppressing star formation and metal enrichment, 
that are also the sources of metal-enriched winds which pre- 
vent Pop III star formation at later times. 

To model this enrichment process in a simple way, we 
assume that a single galaxy launches a wind from the cen- 
ter of each H II region at redshift z ~ 20, and we take it 
that each H II region encloses a comoving volume Vhii = 1 
Mpc 3 h~ 3 at this redshift. While Shin et al. (2008) find that 
the abundance of cosmological H n regions having comoving 
volumes of ^ 0.03 Mpc 3 h~ 3 is ~ 10 times higher than these 
larger H n regions at z ~ 21, overall the larger H II regions 
contain the majority of the reionized volume; thus, we take 
the larger volume as a typical value. We assume that the 
metal-enriched wind propagates outward from the central 
galaxy at a velocity v w j n d, and we compute the comoving 
distance i? W md that it travels, as a function of redshift z: 

rz j. 
#wind(z) = / «wind(l + z')—dz . (6) 

Finally, we estimate the fraction /wind of the H II regions 
that are enriched by the winds as 

/windM ~ Y^gf , (7) 

where we assume that the H n regions each enclose a co- 
moving volume of Vhii = 1 Mpc 3 h~ :i , as explained above. 
To calculate the number density of Pop III star clusters tak- 
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ing into account this metal enrichment by winds, we take 
it that a fraction /wind of haloes that are not self-enriched 
are instead enriched by winds; thus, the number density of 
metal- free star clusters forming in reionized regions nm(z), 
as defined in equation (5), is reduced by this factor. Note 
that for simplicity we have assumed that each individual 
H II region has the same size and launches a wind starting 
at the same redshift; thus the metal-enriched wind volume- 
filling fraction of an individual H II region, given by equation 
(7), is equal to the global fraction in this calculation. 

The resulting number densities are shown in Figure 4. 
The velocity at which metal-enriched winds propagate into 
the IGM at high redshift is not well-constrained, and so we 
plot the number density of Pop III clusters that we obtain 
for several choices, « w ind = 30, 60, and 100 km s _1 . For 
wind velocities of ~ 30 km s _1 , the total number of clusters 
within the DWS area is only slightly lower than in the case 
of no external metal enrichment. For higher velocities the 
expected number of clusters which may be detected drops 
dramatically, with the formation of Pop III clusters being 
prevented altogether below redshifts z ~ 3 and 5, for wind 
velocities of 60 and 100 km s _1 , respectively. To take a fidu- 
cial case, Furlanetto & Loeb (2003) find that the average 
wind velocity from the relatively massive galaxies that we 
are considering may be ~ 60 km s _1 , in which case Pop III 
cluster formation is unlikely to take place at z 3, although 
the abundance of Pop III clusters above this redshift is only 
a factor ^ 2 below the upper limits presented in Section 3.1. 

We emphasize that the simple estimates presented here 
do not take into account processes such as the mixing of 
metals with the primordial gas in the complex cosmological 
density field or the anisotropic expansion of early cosmo- 
logical H II regions. In particular, we note that the mixing 
of metals with the gas contained in pristine haloes may be 
incomplete, leaving the star-forming gas in the centers of 
the haloes pristine, even if the halo as a whole is overtaken 
by a metal-enriched wind (see Cen & Riquelme 2008; see 
also Wyithe & Cen 2007). Furthermore, the metallicity of 
the gas in primordial haloes that is mixed with metals may 
not exceed the threshold necessary for the transition to low- 
mass, Pop II star formation, generally taken to be between 
10~ 6 and 10~ 3 ' 5 Z© (e.g. Bromm & Loeb 2003; Omukai et 
al. 2005; Frebel et al. 2007). In this case, the stellar initial 
mass function (IMF) in these clusters may be very similar 
to the purely metal-free case. However, observations of high 
redshift quasars suggest that much of the enriched IGM at 
z 5 has metallicity Z <; 10~ 4 Z© (Songaila 2001; Pettini et 
al. 2003), such that chemical enrichment by galactic winds 
may indeed effect the transition from a Pop III to a Pop II 
IMF (see also Ferrara et al. 2000; Scannapieco et al. 2003). 
In principle, though, star formation could also be prevented 
entirely in haloes with masses ~ M m i n due to the stripping 
of gas by winds with velocities exceeding the escape velocity 
of the halo (see Scannapieco et al. 2000). 



4 THE DETECTABILITY OF POP III STAR 
CLUSTERS AFTER REIONIZATION 

While we have shown that there may exist a population of 
metal-free haloes in which Pop III star clusters form even af- 
ter reioinization, the detectability of these clusters depends 



on their luminosities and spectral properties. We consider 
the luminosities of the stellar clusters in two emission lines, 
Lya and He II A1640, prominent recombination lines emit- 
ted by ionized hydrogen and doubly ionized helium, respec- 
tively. The ratio of the luminosities of these lines, as well 
as their equivalent widths, are sensitive to the stellar IMF 
and metallicity. In particular, for a given stellar mass, metal- 
free stars are expected to have higher surface temperatures 
than metal-enriched stars, leading to a relatively high ra- 
tio of He n-ionizing photons to H i-ionizing photons being 
emitted. In turn, this leads to relatively strong nebular emis- 
sion in He II A1640 (e.g. Tumlinson & Shull 2000; Bromm et 
al. 2001; Oh et al. 2001; Schaerer 2002, 2003). Detection of 
both Lya and He II A1640 can thus, in principle, facilitate 
the identification of Pop III star clusters (e.g. Tumlinson et 
al. 2001; Nagao et al. 2008; Prescott et al. 2009). More mas- 
sive stars also emit higher numbers of He n-ionizing photons, 
meaning that detection of these lines can also allow to place 
constraints on the stellar IMF. 

While in neutral regions of the IGM Lya is easily scat- 
tered (e.g. Loeb & Rybicki 1999), the star clusters that we 
consider reside in regions which are already reionized, in 
principle yielding the IGM optically thin to Lya, although 
even a small neutral fraction can result in substantial ab- 
sorption (e.g. Gunn & Peterson 1965). In this case, other 
emission lines from ionized hydrogen will still be directly 
detectable, Ha being the strongest of these lines with a lu- 
minosity ~ 10 times lower than that of Lya (e.g. Oster- 
brock & Ferland 2006), and these lines could alternatively 
be used to constrain the metallicity and IMF (see e.g. Oh 
et al. 2001). For our purposes, we will assume that the IGM 
is optically thin to Lya when evaluating the detectability of 
this emission line. 

We consider Pop III star clusters that are formed 
promptly, within the ~ 2 Myr lifetime of the most massive 
individual stars (see also Oh et al. 2001), and we similarly 
evaluate the detectability of these clusters within ~ 2 Myr 
of their formation, as it is during this time that a cluster 
is at its most luminous. To estimate the mass of the stel- 
lar clusters we take it that a fraction /« of the gas that is 
able to collapse into a DM halo is converted into stars; this 
quantity is not well-constrained, but is likely of the order 
of a few percent (e.g. Krumholz & Tan 2007). In turn, we 
assume that only a fraction / co ii ~ 0.1 of the gas associ- 
ated with a DM halo is able to cool and collapse (i.e. that 
the collapsed baryon fraction in the halo is 10 percent of its 
value in the absence of photoheating) , as the star-forming 
haloes under consideration have just grown massive enough 
to allow the photoionized gas to collapse for the first time, 
meaning that the majority of the gas stays photoionized and 
at low densities (e.g. Thoul & Weinberg 1996; Dijkstra et al. 
2004) . We thus estimate the luminosity emitted by a Pop III 
cluster in Lya as 

T in« „ „ -1 ( 1 - /esc, HI \ / /» W/colA 

Ll - ~ 10 crgs { — i — JWJUrJ 

/ Qui \ ( Mmin \ 

X ^1O« S -1M -V 1.1O 8 M J ' W 

where / CS c,hi is the fraction of hydrogen-ionizing photons 
that escape the halo and so do not contribute to the nebular 
emission from the cluster. The number of hydrogen-ionizing 
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Figure 5. The flux in Lyce (top panels) and He II A1640 (bottom panels) from Pop III star clusters. As in Figs. 2 and 3, the solid lines 
correspond to a minimum halo circular velocity of 20 km s — 1 , while the dashed lines correspond to a minimum of 30 km s — 1 ; for each 
of the minimum halo masses, the top (blue) curves corresponds to a top-heavy IMF, while the bottom (red) curves corresponds to a 
Salpeter IMF. The left panels show the flux of clusters assuming a star formation efficiency /* = 10~ 2 , while the right panels correspond 
to the case of /* = 10 — 1 . The other parameters characterizing the clusters are given by our fiducial choices of / esc = and / 00 u = 0.1. 
The black dotted lines in each panel show the flux limits for a 3 a detection in a 100 hr exposure for NIRSpec (ij=100 mode), which will 
be aboard the JWST. Lyce emission will appear in the wavelength range covered by NIRSpec at z > 4, while He II A1640 emission will 
fall in this wavelength at z > 2.7. Detection of Lya would be feasible for all clusters with a top-heavy IMF out to at least z ~ 8, and for 
clusters with a Salpeter IMF and a high star formation efficiency out to at least z ~ 6. However, detection of He II A1640 is only feasible 
for clusters with a top-heavy IMF. The green dotted lines in the bottom panels show the flux limits of the Nagao et al. (2008) survey for 
He II A1640 emitters at z ~ 4; only the most luminous clusters would have been detectable in this survey. 



photons emitted per second per unit mass in stars, Qm, de- 
pends on the IMF, and we consider two possible values: Qhi 
= 10 47 s _1 M Q _1 and 1.6 x 10 48 s _1 M _1 , for a Salpeter 
Pop III IMF and a top-heavy Pop III IMF with a charac- 
teristic stellar mass <; 100 Mq, respectively (Tumlinson & 
Shull 2000; Bromm et al. 2001). As well, we have accounted 
for the fact that only a fraction of ~ 2/3 of the ionizing 
photons that are absorbed are converted into Lya photons 
(e.g. Spitzer 1978). Following the same steps, we have for 
the luminosity in He II A1640 

T A w ln38 -1 / 1 — /csc.HcII \ f f* \ //coll\ 

L 1640 ~ 4X10 ergs ( { j ( — j j 

5 x 10 45 s^M©- 1 J [Wm^ J ' { ' 



where /cscHcii is the escape fraction of He n-ionizing 
photons. For the number of He n-ionizing photons emit- 
ted per second per unit mass in stars, we have QhcII = 
5 x 10 45 s- 1 Mq 1 and 3.8 x 10 47 s _1 M" 1 , again for a 
Salpeter Pop III IMF and a top-heavy Pop III IMF, respec- 
tively (Tumlinson & Shull 2000; Bromm et al. 2001). The 
flux in Lya observed at z = is then given by 

f _ Lj^ytx 

/LyQ = A^Dl 

where Dl{z) is the luminosity distance to redshift z. The 
flux in He n A1640 is computed similarly, by replacing the 
Lya luminosity with the He n A1640 luminosity. 
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Figure 6. The monochromatic Lyo flux of Pop III clusters as would be observed by the NIRCam instrument in the JWST DWS, as 
given by equation (11), for different choices of the IMF, star formation efficiency /*, and minimum halo mass M m ; n , as described in Fig. 
4. The black dotted lines are the 3 a detection limits expected for the DWS, for an exposure time of 2 X 10 s seconds; the different line 
segments correspond to the different wide filters the wavelength range of which contains Lya at redshift z. For a star formation efficiency 
of /* = 10~ 2 (left panel) only clusters with a top-heavy IMF may be detected, while for a high star formation efficiency /* = 10 _1 (right 
panel) even clusters with a Salpeter IMF may be detectable out to z ~ 6. 



In Figure 5 we show the fluxes in Lye* and He II A1640 
from clusters forming in haloes with mass ~ M m i n , for each 
of our choices of this minimum halo mass, and for our two 
choices for the IMF. We also plot the flux limits for 3 a de- 
tection of these emission lines with 100 hours of observation 
by the Near-Infrared Spectrograph (NIRSpec) that will be 
aboard the JWST. NIRSpec will operate at observed wave- 
lengths from 0.6 to 5 fj,m (e.g. Gardner et al. 2006). The 
lower limit of this wavelength range yields NIRSpec sensi- 
tive only to Lye* and He II A1640 at redshifts 2^4 and z i; 
2.7, respectively. 

For the most luminous clusters that we consider, the 
case of a large fraction /* = 10 _1 of the collapsed gas going 
into stars with a top-heavy IMF, both Lya and He II A1640 
would be detectable by NIRSpec out to beyond redshift z 
= 10, as shown in the panels on the right in Fig. 5. Indeed, 
the detection of these emission lines is much more feasi- 
ble for the case of a top-heavy IMF, for all combinations 
of minimum halo mass and star formation efficiency /». In 
particular, because even a 3 a detection of He II A1640 is 
not feasible if clusters have a Salpeter IMF, as shown in the 
bottom panels of Fig. 5, detection of this emission line would 
suggest at least a somewhat top-heavy IMF. We note, how- 
ever, that the escape fractions of both H I- and He n-ionizing 
photons are likely to be higher for the case of a top-heavy 
IMF, owing to the strong hydrodynamic response of the gas 
to photoheating (e.g. Johnson et al. 2009; see also Wise & 
Cen 2009; Razoumov & Sommer-Larsen 2009), with a cor- 
responding decrease in the emitted luminosity. Even so, the 
flux in He II A1640 will likely still, in general, be higher for 
a more top-heavy IMF (Johnson et al. 2009). 

In the bottom panels of Fig. 5 we also plot the approx- 
imate flux limit of the survey for He II A1640 emitters at 
z ~ 4 carried out by Nagao et al. (2008). While this flux 
limit is an order of magnitude higher than that we find for 



NIRSpec, the non-detection of He II A1640 emitters in this 
survey already places some constraint on the abundance of 
Pop III clusters with a top-heavy IMF and a high star for- 
mation efficiency (see Section 3). The JWST DWS will place 
stronger constraints on both the abundance and luminosity 
of Pop III clusters. 

As the DWS will be conducted using the near-infrared 
imaging instrument NIRCam, we would also like to evaluate 
the detectability of the clusters by comparing their expected 
fluxes to the flux limits with which NIRCam will operate for 
the DWS. Given the rest-frame spectral energy distribution 
(SED) of the cluster, the flux that NIRCam would detect 
at z — can be computed. For the expected SED of a clus- 
ter of young Pop III stars the luminosity in Lya is roughly 
an order of magnitude greater than the stellar and nebu- 
lar continuum luminosity in the wavelength range covered 
by a given NIRCam filter (assuming a small escape fraction 
/esc)(e.g. Bromm et al. 2001; Schaerer 2002). For simplicity, 
we thus calculate the flux that would be detected by NIR- 
Cam assuming that it is entirely in Lya. Similar to the cal- 
culation in equation (10), we thus find the monochromatic 
Lya flux to be 

I/LyaALya(l + z)R 

„ 10- 2 nJyf A ^ -Vl±frV^) (11) 
y \10 40 erg s- 1 J V 10 J V 4 / V ' 

where R is the spectral resolution; for the wide filters on the 
NIRCam which operate at the wavelength of Lya, R = 4 §. 
Both the monochromatic Lya flux and the flux limits for the 
planned DWS are shown in Figure 6. While NIRCam has an 

§ http:/ /www. stsci.edu/jwst/science/sensitivity/ 
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exquisite sensitivity which will allow for detection of fluxes 
as low as ~ 1 njy, Pop III clusters may only be detected in 
the DWS if their IMF is top-heavy or if the star formation 
efficiency is relatively high, i.e. of the order of /» ~ 10 _1 . 
Nonetheless, that these clusters would form at lower redshift 
and in more massive DM haloes suggests that their detection 
may be still more feasible than that of Pop III stars in the 
first dwarf galaxies at 2 J; 10 (see e.g. Johnson et al. 2009). 

As discussed in Section 3.2, the fact that these clusters 
would only form within the H 11 regions surrounding the first 
galaxies at 2 J; 20 (which reside in haloes with masses Jj 5 x 
10 7 M ; Trac & Cen 2007) implies that they will be found 
preferentially in the vicinity of the descendants of these first 
galaxies, which the extended Press- Schechter formalism dic- 
tates would typically reside in DM haloes with masses of ~ 
10 11 Mq at 3 ^ 2 < 6. In particular, according to Shin et al. 
(2008) the comoving radius of the largest H 11 regions at z 
^ 20 is ^ 1 Mpc ft -1 . Therefore, a natural prediction of our 
model for Pop III star formation after reionization is that 
Pop III clusters may be found within a comoving distance 
of ~ I Mpc hT 1 , corresponding to an angular separation of 
40 arcsec on the sky, of the galaxies residing in ~ I0 11 
M haloes at 3 2 6. 



5 SUMMARY AND DISCUSSION 

We have shown that the suppression of star formation in DM 
minihaloes by photoheating during the early stages of reion- 
ization can lead to the existence of unenriched DM haloes 
with masses <; 10 8 Mq at redshifts 2 10 which may host 
Pop III star clusters. While the abundance of such haloes is 
sensitively dependent on the degree to which galactic winds 
enrich the IGM with metals, as well as on the minimum 
halo mass for star formation in the reionized regions of the 
Universe, there may be up to thousands of bright Pop III 
star clusters per square degree at redshifts 3 2 <^ 6, corre- 
sponding to ~ 100 clusters within the area of the DWS to be 
carried out by the JWST. This is, however, likely a strong 
upper limit to the abundance, as metal-enriched winds ex- 
pelled by the same galaxies that reionize the primordial gas 
at early times may easily enrich it by these redshifts, depend- 
ing on the speed of the winds and on how quickly the metals 
are mixed into the pristine haloes. If this enrichment is suf- 
ficiently rapid, Pop III star formation in pristine DM haloes 
may be completely prevented after reionization is complete 
at 2 ~ 6. 

Although the stellar mass, and so the luminosity, of 
Pop III clusters formed in reionized regions of the Universe 
is limited by the fraction of the primordial gas that is first 
able to collapse into the DM haloes, we estimate that Lya 
emission from these star clusters could be detected by the 
JWST if the stellar IMF is top-heavy or if the star forma- 
tion efficiency /* is high, i.e. of the order of 10 percent. The 
He 11 A1640 emission line, which is an indicator of stellar 
IMF and metallicity, may be detectable out to at least 2 ~ 
3 for the case of a top-heavy IMF, perhaps allowing to place 
constraints on the Pop III IMF, at least at low redshift. 
Furthermore, in reionized regions Pop III clusters may form 
preferentially within ~ 40 arcsec of more luminous galaxies 
at 3 2 6, the descendants of the first galaxies which 



began reionization at 2 <; 20, perhaps aiding in the identifi- 
cation of these clusters in future surveys. 

Current observations already yield some constraints on 
the abundance of Pop III star clusters at 2 ~ 4 (Dawson et 
al. 2004; Nagao et al. 2008; Wang et al. 2009; see also e.g. 
Fosbury et al. 2003; Shapley et al. 2003; Nagao et al. 2005; 
Jimenez & Haiman 2006; Prescott et al. 2009; Bouwens et al. 
2009). The non-detection of strong He 11 A1640 emitters at 2 
~ 4 reported by Nagao et al. (2008) yields an upper limit on 
the number density of Pop III clusters of nm ~ 2500 Gpc , 
between the upper limits plotted in Fig. 3. However, the 
limiting flux in this survey is only low enough to rule out the 
presence of the most luminous Pop III clusters. Nonetheless, 
this imposes a useful constraint, and suggests that if clusters 
of massive Pop III stars do form after reionization, they 
may in fact form in relatively low-mass haloes in which star 
formation is inefficient. 

We note that the Pop III IMF in star clusters formed 
after reionization is likely different than that of the first stars 
formed in minihaloes at 2 <; 20 (so-called Pop III. I stars), 
due to the onset of turbulence in the more massive haloes 
in which the clusters form (Wise et al. 2007b; Greif et al. 
2008), to the possible enhanced abundances of the coolants 
H2 and deuterium hydride (HD) in the reionized primordial 
gas (e.g. Nagakura & Omukai 2005; Johnson & Bromm 2006; 
Yoshida et al. 2007; McGreer & Bryan 2008), and to the 
lower temperature of the cosmic microwave background (e.g. 
Larson 1998; Tumlinson 2007; Smith et al. 2009). Instead of 
the characteristic stellar mass of ~ 100 Mq expected for 
Pop III. 1 stars formed in the earliest minihaloes, Pop III. 2 
stars with a characteristic mass of the order of 10 Mq may 
instead be detected in star clusters formed from reionized 
primordial gas. Such clusters would have spectral properties 
intermediate between the two IMFs, Salpeter and top-heavy, 
discussed in Section 4. 

An important next step to better estimate the abun- 
dance of Pop III star clusters will be to carry out cos- 
mological simulations which capture the earliest stages of 
the inhomogeneous reionization of the Universe, as well as 
the chemical enrichment from both the first stars formed 
in minihaloes and large-scale galactic winds. We note that 
the cosmological metal-enrichment simulations presented in 
Tornatore et al. (2007) already provide some indication that 
Pop III clusters may form in rare regions after reionization. 
In addition, cosmological radiative transfer simulations may 
be necessary to estimate the masses, and so the luminosities, 
of stellar clusters that form in the DM haloes into which the 
reionized primordial gas is first able collapse in the course of 
hierarchical structure formation. Related to this, magnetic 
fields pervading the reionized gas may have an important 
dynamical effect on the formation of Pop III star clusters at 
low redshifts (see Schleicher et al. 2009a) ; simulations track- 
ing the growth and evolution of magnetic fields in the IGM 
would thus also be of interest. Furthermore, in the present 
work we have only considered the 'average' merger history 
of a potential Pop III star-forming halo, ignoring the range 
of possible merger histories that such haloes may have. This 
assumption would clearly be relaxed within the context of a 
cosmological simulation. 

Although we have accounted for the higher average den- 
sity of regions which are reionized at early times with our 
choice of /reion.mass in equation (5), we have otherwise as- 
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sumed that the haloes in which Pop III star clusters form 
are not biased, i.e. that these haloes are otherwise located 
randomly in space. In principle this could lead to an overes- 
timate of nm(2:), for example, if haloes with mass M m i n (z) 
typically form in underdense regions, while it is only over- 
dense regions which are reionized at early times. However, 
haloes with mass M m i n (z) correspond to J> 1-a fluctuations 
in the cosmological density field at z <; 3 (e.g. Barkana & 
Loeb 2001), typically form at 3 Si 2 S> 4 (see e.g. Boylan- 
Kolchin et al. 2009), and so it does not appear likely that 
they form preferentially in underdense regions at z <; 3, the 
redshifts at which we find Pop III star clusters are mostly 
likely to form in them. Therefore, we expect that accounting 
for the detailed clustering properties of these haloes around 
the descendants of the sources of reionization at z ~ 20 
would not change our results substantially, although the ef- 
fect of bias should be considered in more depth in future 
studies. 



Another implicit assumption throughout this work has 
been that stellar clusters will form in the pristine DM haloes 
which are massive enough for the gas to collapse into them 
after reionization. There is the possibility that this gas could 
instead collapse to form a 10 4 Mq black hole directly, if 
molecular cooling is suppressed (e.g. Bromm & Loeb 2003; 
Begelman et al. 2006; Lodato & Natarajan 2006; Spaans & 
Silk 2006; Regan & Haehnelt 2009; see also Schleicher et al. 
2009b). Shang et al. (2009) find that the LW background at 
z Ss 8 (Ahn et al. 2009) may be high enough to suppress 
the formation of H2 molecules and allow for the formation 
of such a massive single object. However, it is not clear that 
the relatively small fraction / co n of the gas which is able to 
collapse is sufficiently high for this to occur. Also, despite the 
presence of a persistant LW background, the enhanced free 
electron fraction in the collapsing reionized gas will boost 
the formation rate of H2 to some extent, perhaps aiding to 
prevent direct collapse to a black hole. 



Although we have explored the impact on the abun- 
dance of Pop III clusters of different choices for the speed 
of metal enrichment by galactic winds, and for the mini- 
mum DM halo mass for star formation after reionization, 
we have considered only one choice for the early reioniza- 
tion history, based on the simulations of Shin et al. (2008). 
As the abundance of unenriched haloes is dependent on the 
fraction of the Universe reionized at redshifts z J> 20, as 
well as on the typical size of cosmological H n regions at 
these redshifts, our results are in principle quite sensitive to 
this choice. Observational constraints on the early stages of 
reionization, for example from observations of 21-cm emis- 
sion by instruments such as the Square Kilometer Array (e.g. 
Carilli et al. 2004), may thus provide indirect constraints 
on the abundance of unenriched haloes at lower redshifts. 
In the upcoming decade, next generation telescopes such as 
the JWST and the E-ELT, however, may first provide direct 
constraints on the amount of Pop III star formation taking 
place in the early Universe, with the exciting possibility of 
detecting the emission from individual Pop III star clusters, 
perhaps even well after the epoch of reionization. 
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